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I. INTRODUCTION 

Due to the high fo-quark mass, weak decays of baryons 
containing a b quark are a good testing ground of some 
approximations in quantum chromodynamics (QCD) cal- 
culations, such as heavy-quark effective theory (HQET) 
[y. Alternatively, when using such calculations, the A° 
may provide a determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) couplings with systematic uncertainties 
different from the determinations from the decays of B 
mesons Q. While the B mesons are well studied, less is 
known about the A b baryon. Only nine decay modes of 
the A b have been observed so far, with the sum of their 
measured branching fractions of the order of only 0.1 and 
with large uncertainties on the measurements [3j. While 
theoretical predictions are available for the A° — > A+7r~ 
branching fraction no prediction is currently avail- 
able for the A b — > A+7r~7r + 7r~ decay mode. LHCb re- 
cently reported the measurement of the ratio of branch- 
ing fractions B(A° b -> A+7r-7r+7r-)/2?(Ag -> A+tt") 
= 1.43 ± 0.16(stat) ± 0.13(syst) [|. 

This paper reports a study of the A b — > A+7r _ 7r+7r _ 
decay mode and is especially distinguished by the high 
yields and high precision measurement of the A b — > 
A+7r~7r + 7r~ resonant contributions, the following decay 
modes: 



A° b ^ A c (2595)+tt-, 
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A° -> A c (2625)+tt-, 
Ag -> S c (2455) ++ 7r-7r-, 
A° -> E c (2455)°7r+7r-. 

We measure the branching fraction of each resonant 
decay mode relative to the A" — > A+7r~ decay 
mode, and the ratio of branching fractions B(A° — > 
A+7r _ 7r + 7r _ )/B(A" — > A+7r~). The measurement is per- 
formed using a sample of pp collisions corresponding to 
2.4 fb _1 integrated luminosity collected by CDF II be- 
tween February 2002 and May 2007. We reconstruct Ag 
decays from particles whose trajectory projections in the 
plane transverse to the beamline do not intersect the 
beamline (displaced tracks). The signal yields of inter- 
est are extracted by fitting mass differences to minimize 
the effect of systematic uncertainties. As a crosscheck, 
we repeat the analysis on the reference decay modes 
B" -> D-7t + t:-tt + and B° -> D~Tr + . 

The structure of the paper is as follows. Section |TT] 
describes the detector systems relevant to this analysis. 
Event selection and A° — > A+7r~7r + 7r~ and A° — > A+7r~ 
candidate reconstruction are described in Sec. IIIII In 
Sec.QV]we present the signal yields. In Sec.|V]we describe 
the evaluation of the detector acceptance and the rela- 
tive branching fraction measurements, while in Sec. IVII 
the systematic uncertainties arc discussed. Final results 
are reported in Sec. IVII1 



II. THE CDF II DETECTOR AND TRIGGER 

The CDF II detector is a multipurpose magnetic spec- 
trometer surrounded by calorimeters and muon detec- 
tors. The components relevant to this analysis are briefly 
described here. A more detailed description can be found 
elsewhere Q. A silicon microstrip detector (SVX and 
ISL) and a cylindrical drift chamber (COT) Q im- 
mersed in a 1.4 T solenoidal magnetic field allow the re- 
construction of charged particle trajectories in the pseu- 
dorapidity @ range \r)\ < 1.0 [T^j. The SVX detector 
consists of microstrip sensors arranged in six cylindrical 
shells around the beamline with radii between 1.5 and 
10.6 cm, and with a total z coverage of 90 cm. The first 
SVX layer, also referred to as the LOO detector, is made 
of single-sided sensors mounted on the beryllium beam 
pipe. The remaining five SVX layers are made of double- 
sided sensors and divided into three contiguous five-layer 
sections along the beam direction z. The two additional 
silicon layers of the ISL help to link tracks in the COT to 
hits in the SVX. The COT has 96 measurement layers be- 
tween 40 and 137 cm in radius, organized into alternating 
axial and ±2° stereo superlayers. The charged particle 
transverse momentum resolution is o Pt /pt — 0.07% pr 
(GeV/c), and the resolution on the transverse distance of 
closest approach of the particle trajectory to the beam- 
line (impact parameter, do) is ~40 /im, including a «30 
/im contribution from the beamline. 



Candidate events for this analysis are selected by a 
three-level on-line event selection system (trigger). At 
level 1, charged particles are reconstructed in the COT 
axial superlayers by a hardware processor, the Extremely 
Fast Tracker (XFT) [ll|. Two charged particles are re- 
quired with transverse momenta pt > 2 GeV/c. At level 
2, the Silicon Vertex Trigger (SVT) [U[ associates SVX 
r — (ft position measurements with XFT tracks. This pro- 
vides a precise measurement of the track impact param- 
eter do- We select 6-hadron candidates by requiring two 
SVT tracks with 120 /im < d < 1000 /im. To reduce 
background from light-quark jet pairs, the two trigger 
tracks are required to have an opening angle in the trans- 
verse plane 2° < A<fi < 90°. The tracks must also satisfy 
the requirement Lt > 200 /im, where Lt is defined as 
the distance in the transverse plane from the beam line 
to the two-track intersection point, projected onto the 
two-track momentum vector. The level 1 and 2 trigger 
requirements are then confirmed at trigger level 3, where 
the event is fully reconstructed. 



III. EVENT RECONSTRUCTION 

The search for A" A+7r~7r+7r~ and Ag A+7r~ 
candidates begins with the reconstruction of the A+ using 
the three-body decay A+ — s- pK~n + [l3j ]. Three tracks, 
assumed to be a kaon, a proton, and a pion, with a total 
charge of +1, are fit to a common vertex. No particle 
identification is used in this analysis. All particle hy- 
potheses consistent with the candidate decay chain are 
considered. Additional selection criteria (cuts) are ap- 
plied on fit probability (P(x 2 (A+)) > 10 -5 ), transverse 
momentum (py(A+) > 4.0 GeV/c), and transverse de- 
cay length relative to the beamline (Lt(A+) > 200 /im). 
We also require prip) > Pt{^ + ), to suppress random- 
track combinatorial background. The reconstructed A+ 
mass (m(A+)) distribution is comparable to the one re- 
ported in Ref. [l4| . The reconstructed A+ mass is re- 
quired to be close to the known A+ mass (2.240 - 2.330 
GeV/c 2 ) 1. Since mass differences arc used to search for 
the resonances, no mass constraint is applied in the A+ 
reconstruction. The A" A+7r~7r+7r _ (Ag A+7r~) 
candidate is reconstructed by performing a fit to a com- 
mon vertex of the reconstructed A+ and three (one) ad- 
ditional tracks, assumed to be pions, with pt > 0.4 
GeV/c, and a total charge of -1. For all the possible 
track pairs out of the six (four) tracks that form the A° 
candidate, we require the difference between the z coor- 
dinate of the points of closest approach of the two tracks 
to the beam to be less than 5 cm. Additional cuts on the 
A° candidate fit probability (P(x 2 (A°)) > 10~ 4 ), trans- 
verse momentum (p^(A°) > 6.0 GeV/c), transverse de- 
cay length relative to the beamline (Lt(A°) > 200 /im), 
and A+ transverse decay length relative to the beam- 
line (L T (A+) > 200 /im) and to the Ag vertex (L T (A+ 
from Aj|) > —200 fim) are applied. We also require that 
the transverse momentum of the pion produced in the 
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A+ decay is larger than the transverse momentum of the 
same-charge pion produced in the Ag decay, which con- 
siderably reduces the combinatorial background due to 
the larger boost of the pion produced in the A+ decay. 
To improve the purity of the Ag — > A+7r~7r + 7r~ signal, 
we optimize the analysis cuts to maximize the signal sig- 
nificance S/y/S + B. The number of Ag — > A+7r _ 7r + 7r~ 
candidates S and the number of background events B are 
estimated in data by performing a fit of the m(Ag) dis- 
tribution. This procedure determines the final selection 
criteria: p T (A° b ) > 9.0 GeV/c, L T (A° b ) / a Lt{K) > 16, 
do(Ag) < 70 /im, and Ai?(7r"7r + 7r") < 1.2, where d (Ag) 
is the impact parameter of the reconstructed Ag candi- 
date relative to the beamline and Ai?(7r~7r + 7r~) is the 
maximum yj Arj 2 + A</> 2 distance between the two pions 
in each of the three possible pairs of pions. We veri- 
fied that by splitting the data sample in two independent 
samples, the optimization procedure yields the same final 
selection criteria when applied separately to the two sam- 
ples, and that the Ag — > A+7r~7r + 7r~ yield is evenly dis- 
tributed. This ensures that our optimization procedure 
does not introduce a bias on the branching fraction mea- 
surement. To reduce possible systematic effects in the 
estimate of the reconstruction efficiency due to Monte 
Carlo simulation model inaccuracy, the same selection 
cuts optimized for A° — > A+7r - 7r + 7r~ are also applied to 
the selection of the Ag 

Ai?(7T~7r + 7T~) CUt. 



the Ag -> A+7T 



A+7r signal, except for the 



DETERMINATION OF THE SIGNAL 
YIELDS 



IV. 



Figure Q [](a)| shows the distribution of the difference be- 
tween the reconstructed Ag and A+ masses, m(Ag) — 
m(A+), of the selected Ag — > A+7r~7r + 7r~ candidates 
with the fit projection overlaid. A significant signal of 
Ag — > A+7r _ 7r + 7r _ is visible centered approximately at 
3.330 GcV/c 2 . Backgrounds include misreconstructed 
multibody 6-hadron decays (physics background) and 
random combinations of charged particles that acci- 
dentally meet the selection requirements (combinatorial 
background). We use an unbinned extended maximum- 
likelihood fit to estimate the Ag — > A+7r~7r + 7r~ signal 
yield. The signal peak is modeled with a Gaussian, with 
mean and width left floating in the fit. The combinato- 
rial background is modeled with an exponential function 



of m(Ag) — m(A^) with floating slope and normaliz£ 
tion. The distribution of the main physics backgrounds, 



due to the B9. 



D 



(*)- 



✓ , * /i 7r 7r + decay modes, are de- 
rived from simulation and included in the fit with fixed 
shape and floating normalization. The Ag — > A+7r~7r + 7r~ 
yield estimated by the fit of the data is 1087±101 can- 
didates, the world's largest sample currently available of 
this decay mode. Figure QP(b)| shows the Ag mass dis- 
tribution of the selected Ag — >• A+7r _ candidates. The 
A? mass distribution is described by several components: 



Gaussian signal, a combinatorial back- 
ground, reconstructed B mesons that pass the A+7r~ se- 
lection criteria, partially reconstructed A® decays (e.g. 
Ag — > A£l~Di), and fully reconstructed Ag decays other 
than A+7r~ (e.g. A® — > A^K~). Also in this case the dis- 
tributions of physics backgrounds are derived from sim- 
ulation and included in the fit with fixed shapes and 
floating normalization, as detailed in Ref. 15|. The 
A" — > A+7T" yield estimated by the fit of the data is 
3052±78 candidates. 

In the reconstructed Ag — > A+7r~7r + 7r~ sam- 
ple we searched for the resonant decay modes: 
A" -> A c (2595)+tt-, Ag -> A c (2625)+tt-, Ag -> 
S c (2455) + +7r~7r-, and Ag -> £ c (2455)°7r+7r-. The 
available energy transferred to the decay products in the 
decays of the charmed baryons (A c (2595) + , A c (2625) + , 
£ c (2455) ++ , and £ c (2455)°) into A+ is small. There- 
fore the differences of the reconstructed masses m(A* + ) — 
m(A+), m(£ c (2455)++) - m(A+), and m(£ c (2455)°) - 
m(A+) are determined with better resolution than the 
masses of the charmed baryons, since the mass resolution 
of the A+ signal and most of the mass systematic uncer- 
tainties cancel in the difference. Figure [ ^a)| shows the 
m(A*+)-m(A+) distribution, for Ag -> A+7r~7r+7r _ can- 
didates with mass in a ±3er range (±57 MeV/c 2 ) around 
the Ag mass. The A c (2595)+ and A c (2625)+ signals are 
clearly visible. Although there are two possible A* + can- 
didates for each Ag — > A+ir~ir + ir~ decay, only the can- 
didate made with the tt~ with lower pt has a value of 
ra(A* + ) — m(A+) in the mass region where the A c (2595) + 
and A c (2625) + signals are expected. The A c (2595) + and 
A c (2625) + signal yields are estimated with an unbinned 
extended maximum-likelihood fit. The A c (2595) + and 
A c (2625) + signals are modeled with two non-relativistic 
Breit-Wigner functions convolved with the same Gaus- 
sian resolution function, since the mass difference be- 
tween the two resonances is tiny. The background is 
modeled by a linear function. The A c (2595) + natural 
width is mass dependent to take into account the thresh- 
old effects, as reported in Ref. 

[3, the A c (2625)+ natural 
width and the width of the Gaussian resolution function 
are free parameters of the fit. Table |T] reports the esti- 
mated signal yields and significances, evaluated by means 
of the likelihood ratio test, LR = L/Lb c k, where L and 
Lbck are the likelihood of the signal and no signal hy- 
potheses, respectively fil| . 

Figures [|fb71 and ffic)] show the m(Ag) - m(A+) distri- 
bution restricted to candidates with m(A* + ) — m(A+) < 
0.325 GeV/c 2 and 0.325< m(A*+) - m(A+) < 0.360 
GeV/c 2 , respectively, i.e. compatible with the A c (2595) + 
and A c (2625) + expected signals. Each signal is mod- 
eled with a Gaussian function, with floating mean and 
width. The combinatorial background is modeled with 
an exponential function with floating slope and normal- 
ization, and the physics background, which is mainly 
due to semileptonic Ag —> A+ir~ir + l~Vi decays, is de- 
rived from simulation and introduced in the fit with fixed 
shape and floating normalization. We verified that the 
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FIG. 1: The reconstructed invariant mass spectra after ap- 
plying all the selection criteria: (a) the mass difference 
m(Ab) — m(Ac) distribution of the A° — > A^7r~7r + 7r _ can- 
didates; (b) the m(A°) distribution of the A° — ¥ A^7r - can- 
didates. 



m(A*+) - m(A+) > 0.380 GeV/c 2 . In Fig. E gg] 
and E fb)1 the resulting m(E c (2455)++) - m(A+) and 
to(E c (2455)°) — m(A+) distributions arc shown. Promi- 
nent S c (2455) ++ and S c (2455)° signals are visible. 
While there is only one S c (2455) ++ candidate for each 
A" — > A+7r _ 7r + 7r~ decay, two S c (2455)° candidates are 
possible. Also in this case, only the candidate made with 
the ir~ with lower pt is in the E c (2455)° mass region. 
The S c (2455) ++ and E c (2455)° signals are modeled with 
non-relativistic Breit-Wigner functions convolved with a 
Gaussian resolution function, with the addition of an 
empirical background [13, GJ]. The E c (2455)++ and 
E c (2455)° natural widths are Gaussian constrained to the 
world average values @ , while the width of the Gaussian 
resolution function is determined to be 1 MeV/c 2 from 
larger statistics samples of E c (2455) ++ and S c (2455)° 
in the A° lower mass region and is fixed in the fit. 
The effect of this approximation is taken into account 
in the systematic uncertainties. The estimated A° — > 
£ c (2455) ++ 7r-7r- and Ag -> E c (2455) 7r+7r~ yields and 
significances arc reported in Tabic |IJ 

In Fig. EfcJ] and Egd)] the m(Ag) - m(A+) distri- 
butions are shown restricted to candidates with 0.160 
< m(£ c (2455) ++ <°) - m(A+) < 0.176 GeV/c 2 , where 
the E c (2455) ++ and S c (2455)° signals are contained. 
The A° signal is modeled with a Gaussian distribution, 
with floating mean and width, while the combinato- 
rial background is an exponential function with floating 
slope and normalization. We verified that the A° — s- 
£ c (2455) ++ 7r-7r- and A° b -> S c (2455)°7r+7r- yields es- 
timated by fitting the m(A^) — m(A+) distributions are 
compatible with the yields reported in Table U with lower 
significance. The fitted masses and widths of the four res- 
onances are in agreement with the world averages Q and 
the recent CDF II measurements UM. 



TABLE I: Yields and significances of the A° — > A+7T ty + -k 
decay modes. The quoted uncertainty is statistical only. 



AjJ decay mode 



— > Ac"K~ir + ir~ 

— > A+7T~7r + 7T~ 

~7r~ — !> A^7r~7r + 7r~ 
E c (2455)°7r + 7r- -> A+tt^tt+tt" 
A^7r~7r + 7T~ (other) 



Yield Significance(cr) 



A C (2595) H 7T 
A c (2625)+tT 
£ c (2455) ++ - 



46.0 ±8.2 
135 ± 15 
110 ±19 
36 ±11 

790 ± 100 



6.2 

>8 
6.6 
3.4 

>8 



A° -> A c (2595)+tt- and A° -> A c (2625)+tt- yields es- 
timated by fitting the m(A®) — m(A+) distributions are 
compatible with the yields reported in Table U with lower 
significance. 

S c (2455) ++ 7r-7r" and A° -> 
signals, the contributions due to the 
A c (2595)+7r" and A° -> A c (2625) + 7r~ decay 
modes are removed by applying the veto requirement 



To extract the A" 
£ c (2455)°7r+7r 



The residual A° signal (named Ag — > 
A+7r~7r + 7r~ (other)) is selected by applying the 
cuts m(A* c +) - m(A+) > 0.380 GcV/c 2 and 
to(S c (2455) ++ '°) - m(A+) > 0.190 GcV/c 2 to re- 
move the contribution due to the resonant decay modes 
(Fig. 0]). This residual A° signal is likely due to a com- 
bination of the A£ -> A+ai(1260) _ , Ag -> A+p°7r" with 
non-resonant p°7r~ (i.e. not produced by a ai(1260)~ 



decay), and non-resonant Ag — > A+7r~ 



decay 



modes, in unknown proportions. A fit is performed with 
a Gaussian function, with floating mean and width to 



7 



model the signal, an exponential function with floating 
slope and normalization to model the combinatorial 
background, and a physics background due to the 
B® s ^ — > D^*) 7r + 7r~7r + decay modes, derived from 
simulation and included in the fit with fixed shape and 
floating normalization. The resulting yield is 790±100 
candidates (Table HJ. The unknown composition of the 
A b — > A+7r~7r + 7r~ (other) sample is taken into account 
as a source of systematic uncertainty. 

V. MEASUREMENT OF THE RATIO OF 
BRANCHING FRACTIONS 

S(A° -> A+7Y-TT + TV-)/B(A° b ->■ A+7T") 

We measure the following ratio of branching fractions: 
B(A°-> A+O ~ 

_ iV(Ag -> I -> A+TT-TT+TT-) e A^A+7r- 

where N are the measured signal yields reported in Ta- 
ble U and the sum on the intermediate V states in- 
cludes A C (2595)+7T-, A C (2625)+7T~, £ c (2455)++tt-tt-, 
S c (2455)°7r+7r-, and A+tt-tt+tt- (other). In the last 
state, we assume equal proportions of the three decay 
modes A° b -> A+Oi(1260) _ , A" b -> A+^tt", and non- 
resonant A° — > A+tt-tt+tt-. To convert event yields 
into relative branching fractions, we apply the correc- 
tions £A°->-A + 7r- / ei ^ or ^ ne var i° us trigger and offline se- 
lection efficiencies of the decay modes A b — > A+tt - and 
A b — > i — > A+7T _ 7r + 7r _ . All corrections are determined 
from the detailed detector simulation. The bgenera- 
TOR program produces samples of specific B hadron de- 
cays according to measured pt and rapidity spectra fl9| . 
Decays of b and c hadrons and their daughters are simu- 
lated using the evtgen package [2(| • The geometry and 
response of the detector components are simulated with 
the GEANT software package [2l[ and simulated events 
are processed with a full simulation of the CDF II de- 
tector and trigger. The resulting estimated corrections 
e A°->AU-/ e i are 4J0 ± °T0, 4.66 ± 0.10, 5.28 ± 0.11, 
and 18.49 ± 0.66, respectively, for the A c (2595) + 7r~, 
A c (2625)+tt-, E c (2455) ++ 7r-7r-, and £ C (2455)°TT+Tr- 
decay modes. For the A+7r~7r + 7r~ (other) decay mode 
a correction factor equal to 9.16 ± 0.14 is obtained by 
averaging the relative efficiencies of the three interme- 
diate states Ag -> A+oi(1260)~, Ag -> A+p°7r", and 
non-resonant A" — > A+ ir~ tt+ tt~ . 

With a similar method, we also measure the ratios of the 
branching fractions of the intermediate resonances con- 
tributing to Ag -> A+7r _ 7r + 7r~, 

BjA" -» j -» A+TT-TT+TT-) _ 
S(Ag -> A+TT-TT+TT-) ~ 



N(A° -> j -> A+TT-TT+TT") 
Ei N(A° i A+TT-TT+TT-)| ' 



VI. SYSTEMATIC UNCERTAINTIES 

The dominant sources of systematic uncertainty are 
the unknown relative fractions of A b —> A+ai(1260)-, 
A b —> A+p°tt _ , and non-resonant A b —> A+tt-tt+tt- , 
which affect the A b — > A+tt^tt+tt - (other) decay mode 
efficiency, and the unknown A b production and A+ 
decay polarizations, which affect the estimate of all 
the £i and e^-^A+Tr- efficiencies. The correction 

e A°->A+,r-/ e A°->A+,r-7r+7r- (other) has an average value of 
9.16 and varies between a minimum of 7.4 and a maxi- 
mum of 11.6, obtained in the extreme cases in which the 
A b — > A+Tr~Tr+Tr~ (other) sample is assumed to be en- 
tirely composed of A b — > A+ai(1260)~ or non-resonant 
A" — > A+7r~Tr+7r~, respectively. The dependence of 
B(A° b -> A+tt-tt+tt- )/B(A° b -> A+TT-) on the frac- 
tion of A° b -> A+ai(1260)- and A° -> A+p -K~ in the 
A° — > A+tt-tt+tt- (other) sample is shown in Fig. [5] The 
difference between the values computed with the average 
and the minimum (maximum) efficiency correction, re- 
spectively, is taken as an estimate of the lower (upper) 
associated systematic uncertainty. 

The unpolarizcd A b and A+ simulation samples are 
used to obtain the central values of the efficiency correc- 
tions. For the study of the systematic uncertainties, an- 
gular distributions in simulation are reweighted accord- 
ing to all possible combinations of the A b production 
polarization states along the normal to the production 
plane, with the A+ polarization states. The A b polar- 
ization and the A+ polarizations arc both taken to vary 
independently in the range ±1. We assume the extreme 
scenarios where both the A b and A+ baryons are 100% 
polarized and we recompute the efficiency corrections as- 
suming the four possible A b and A+ polarization combi- 
nations. The difference in the efficiency corrections be- 
tween the simulation with reweighted angular distribu- 
tions and the simulation with unpolarizcd A b and A+ 
is used to determine the associated systematic uncer- 
tainty. These two sources of systematic uncertainty ac- 
count for approximately 98% of the total systematic un- 
certainty on the measurement of the relative branching 
fraction B(A° b -> A+tt-tt+tt- )/B(Ag -> A+tt"). Other 
systematic errors stem from the uncertainties on the 
A b —> A+tt" background shapes; on the Cabibbo sup- 
pressed decay mode contributions, which affect the esti- 
mate of the signal yields; on the Monte Carlo simulation 
of the signal decay modes (limited sample statistics, trig- 
ger emulation, and Ag production transverse momentum 
distribution), which affect the estimate of the efficiency 
corrections. The contributions due to the uncertainties 
on the £++ and £° signal and background shapes, the 
A+ and A*+ branching fractions, and the A b and A+ 
lifetimes are negligible. 
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FIG. 2: The A° -s- A c (2595) + tt" and A° -)■ A c (2625) + tt~ signals: (a) m(A*+) - m(A+) distribution for candidates in a ±3cr 
range (±57 MeV/c 2 ) around the A^ mass; (b) m(A°) — m(A+) distribution restricted to candidates in the region m{A^ + ) — 
m(A+) < 0.325 GeV/c 2 ; (c) m(A°) — m(Aj) distribution restricted to candidates in the region 0.325 < m(A*+)-m(A+) < 0.360 
GeV/c 2 . 



As a cross-check of the analysis, we also mea- 
sure the relative branching fraction B(B° — > 
D~TT + Tr~ir + )/B(B —> D~ir + ), using the same data 
sample and vertex reconstruction procedure developed 
for the A° analysis. We apply the same optimized cuts 
to the B° candidates, with the additional request to have 
a D~ candidate with mass within ±22 MeV/c 2 of the 
known mass of D~ We estimate B° — > D~7r + 7r~7T + 
and B° -> D~tt + yields of 431 ± 32 and 1352 ± 44 
candidates, respectively. Our measurement B(B° — > 



D-TT + TT-TT + )/B(B° -> D-7T+) = 3.06 ± 0.25(stat) is 
in good agreement with the value calculated from the 
measured absolute branching fractions of the B° decay 
modes reported in Ref. 
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TABLE II: Measured branching fractions relative to the A° — > A+7r~ decay mode (second column). Absolute branching 
fractions (third column) are derived by normalizing to the known value B(Al A^7r~) = (8.8 ± 3.2) x 10~ 3 [22j|. The first 
quoted uncertainty is statistical, the second is systematic, and the third is due to the uncertainty on the A° — > A^tt~ branching 
fraction. 



AjJ decay mode Relative B to Aj| Aj tt~ Absolute B(10~ 3 ) 
B(A b ^ A c (2595) + 7r~) ■ B(A C (2595)+ ->■ A+tt+tt") (7.1 ± 1.3 ±0.6) • 10~ 2 0.62 ± 0.11 ± 0.05 ± 0.23 
B(A° -)■ A C (2625) + 7T-) ■ B(A C (2625) + -> A+tt+tt") (20.6 ± 2.4±^|) ■ 10" 2 1.81 ± 0.21±g;i| ± 0.66 
6(A£ ->■ E c (2455) ++ 7r-7r") • S(E C (2455)++ ->■ A+tt+) (19.0 ± 3.3 ±1.1) • 10" 2 1.67 ± 0.29 ± 0.10 ± 0.61 
B(A° b -> E c (2455)°7r + 7r-) • S(E C (2455)° -J- A+tT) (21.5 ± 6.5±^) • 10" 2 1.89 ± 0.57±^° ± 0.69 
5(A£ -)■ A+7r~7r+7r^ (other)) 2.36 ± 0-32±g':f| 20.8 ± 2.8i|; ( 7' ± 7 - 6 

g(Ag A+tt-tt+tt-) 3.04 ± 0.33lg;^ 26.8 ± 2.9±|;| ± 9.7 



VII. RESULTS 

We measure the relative branching ratio of A° — ^ 
A+7r _ 7r + 7r~ to A° — > A+7r~ decays to be 

S|z^p.3.04 ± 0.33 (s ta t , iS: ™, syst) . 

The relative branching fractions of the intermediate 
states contributing to A° — > A+7r~7r + 7r~ with respect 
to A° — >■ A+7T - are reported in Table [TTJ The absolute 
branching fractions are derived by normalizing to the 
known value £(A° -> A+tt") = (8.8 ± 3.2) x 10" 3 
To compare our result with the recent LHCb measure- 
ment Q of 1.43 ± 0.16(stat) ± 0.13(syst), we assume 
the composition of the admixture to be 2/3 A" — > 
A+ai(1260) - and 1/3 Ag -> A+^tt", and use the overall 
A° — > A+7r~7r + 7r~ yield and a global efficiency correc- 
tion to compute B{A% -> A+7r~7r+7r~)/,B(Ag -> A+7r~), 
as in the LHCb analysis. This results in a value of 
2.55 ± 0.25(stat)j^Q'27(syst), which is inconsistent with 
the LHCb result at the level of 2.6 Gaussian standard 
deviations. 

We also measure the relative branching fractions of 
the intermediate resonances contributing to the Ag — > 
A+7r~7r + 7r~ decay (Table Ullf . These results are of com- 
parable or higher precision than existing measurements. 



VIII. CONCLUSION 

In summary, we reconstruct the A° — > A+7r~7r + 7r~ 
decay mode and the A° — > A c (2595) + 7r~, 
A° -> A C (2625)+7T-, A° -> E c (2455)++7r-7r-, 



and Ag — > E c (2455)°7r + 7r~ resonant decay modes 
in CDF II data corresponding to 2.4 fb _1 of in- 
tegrated luminosity. We measure the branching 
fraction of the resonant decay modes relative to 
the A° — > A+7r~ branching fraction. We also 
measure B(A° b -> A+7r'7r+7r-)/6(Ag -> A+tt") = 
3.04 ± 0.33(stat)tH2(syst). Using the known value of 
B(A° b -> A+7T-) we find £(A° -> A+7i-7r + 7r-) = 
(26.8 ± 2.9(stat)t^(syst) ± 9.7(norm)) x 1CT 3 , where 
the third quoted uncertainty arises from the A° — > A+7r~ 
normalization uncertainty. 
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FIG. 4: The Ag ->■ A+ir~ir + ir~ (other) signal after vetoing the 
resonant decay modes: m(Aj) — ra(A+) distribution. 
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FIG. 5: B{A° b -> A+n-TY+n-)/B(A° b A+tt") (color 
scale) as a function of the assumed fractions of A° — > 
A^a^ and A° — > A+p°7r _ in the composition of the 
A° — > A^7r _ 7r + 7r _ (other) sample. The central value of 
the ratio is overlaid in each bin. The fraction of non- 
resonant AjJ — > A^7r~7r + 7r _ is equal to 1 — / (A° — > A+a^) — 

/ (A° —y A^p ^ - ) . The cross represents the composition 
chosen for the present measurement assuming equal propor- 
tions of A" — s> A+a^, A[J — > A+p°7r~ and non-resonant 
A° -)• A+7r~7r + 7r~. 



